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ABSTRACT 

We report on a 90 ks Suzaku observation of the radio-loud quasar 4C +74.26. The 
source was observed in its highest flux state to date, and we find that it brightened 
by about 20 per cent during the observation. We see evidence of spectral hardening as 
the count rate increases and also find that the rms variability increases with energy up 
to about 4 keV. We clearly detect a broadened Fe line but conclude that it does not 
require any emission from inside about 50 rg, although a much smaller inner radius 
cannot be ruled out. The large inner radius of our best fit implies that the inner disc 
is either missing or not strongly illuminated. We suggest that the latter scenario may 
occur if the power-law source is located high above the disc, or if the emission is 
beamed away from the disc. 

Key words: galaxies: active - galaxies: individual: 4C -f 74.26 ~ X-rays: galaxies. 
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1 INTRODUCTION 

Only a minority of active galactic nuclei (AGN) harbour 
relativistic radio-jets. The mechanism that drives the jet 
production in these radio-loud objects is unknown and con- 
stitutes an important issue in the study of accreting black 
holes. As both jets and X-ray emission originate close to 
the central black hole, it is clear that insight into this prob- 
lem can be gained by comparing the X-ray properties of 
radio-loud and radio-quiet objects (see Ballantyne 2007 for 
a recent review). 

Such studies of radio-loud and radio-quiet AGN, as 
well as Galactic black holes in their different states, have 
shown that radio-loud objects generally have weaker reflec- 
tion features (as measured by the Fe Ka line and Compton 
hump) and narrower Fe lines than their radio-quiet coun- 
terparts (see e.g Wozniak et al. 1998; Eracleous, Sambruna 
& Mushotzky 2000; Grandi, Malaguti & Fiocchi 2006 for 
AGN and Fender, Belloni & Gallo 2004 for Galactic black 
holes). Explanations for these findings include highly ion- 
ized accretion discs (Ballantyne, Ross & Fabian 2002) and 
dilution by jet emission (e.g. Page et al. 2005), but the usu- 
ally favoured interpretation is that the thin accretion disc in 
radio-loud objects is truncated at large radii and replaced by 
an advection-dominated accretion flow (ADAF). The latter 
explanation is especially appealing in terms of jet produc- 
tion, as ADAFs are thought to naturally lead to bipolar 
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outflows (Narayan & Yi 1994, 1995). While ADAFs are ex- 
pected to exist within objects with very low accretion rates 
(M/MEdd < 0.01, Rees et al. 1982), such as low-luminosity 
AGN and Galactic black holes in their low/hard states, they 
are unlikely to be present in radio-loud quasars, which gen- 
erally have higher accretion rates. 

Over the past few years there have been several re- 
ports of radio-loud objects, most of them quasars, where 
the presence of very broad Fe lines suggest that the accre- 
tion disc extends close or all the way in to the innermost sta- 
ble orbit. Examples include the radio-loud AGN 4C +74.26 
(Ballantyne & Fabian 2005), 3C 273 (Tiirler et al. 2006), 
PG 1425+267 (Miniutti & Fabian 2006), 3C 109 (Allen et 
al. 1997; Miniutti et al. 2006) and 3C 120 (Kataoka et al. 
2007), as well as the Galactic black hole GX 339-4 in its 
low/hard state (Miller et al. 2006). These findings strongly 
indicate that a truncated thin disc is not a necessary condi- 
tion for jet formation, and in particular not appropriate for 
objects with relatively high accretion rates. In this paper we 
will present the analysis of a new Suzaku observation of one 
of these objects, 4C +74.26. 

4C +74.26 is a broad line radio galaxy (BLRG) located 
Sit z — 0.104 (Riley et al. 1989). Its radio luminosity is on the 
boarder of FRII and FRI but its morphology is clearly that 
of an FRII object (Riley et al. 1989). A one-sided jet which 
is at least 4 kpc long has been observed with the VLA (Riley 
& Warner 1990), and on parsec scales with VLBI (Pearson et 
al. 1992) . Based on the flux-limit for a counter jet, Pearson et 
al. (1992) find that the inclination of the source axis must be 
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less than 49° to the Hne of sight. The bolometric luminosity 
of the source is around 2 x W^^ erg s~^, making it a low- 
luminosity quasar (Woo & Urry 2002). The same authors 
also estimate a black hole mass of 4 X lO'' Mq, which means 
that the source is operating at an Eddington ratio of ~ 0.04. 

4C +74.26 was first detected in the X-rays in the 
ROSAT All-Sky Survey, but the first observation that al- 
lowed a detailed spectral analysis was carried out with 
ASCA in 1996 (Brinkmann et al. 1998). Brinkmann et al. 
(1998) detected an Fe line with EW ^ 100 eV but could not 
determine if the line was broadened. The continuum could 
be well fitted with a model consisting of a power law and cold 
reflection (with a large reflection fraction of ^ 6), modified 
by absorption at low energies. The ASCA data were later 
re-analysed by Sambruna, Eracleous & Mushotzky (1999), 
who favoured a double power-law model for the continuum, 
in which the hard power law was attributed to emission from 
a jet. For the double power-law model the Fe line was found 
to be broad (cr ~ 0.6 keV) and the equivalent width higher 
(~ 200 eV). 

The spectrum of a subsequent _BepposSAX observation 
could be described with a similar model as the one presented 
by Brinkmann et al. (1998), but with a significantly smaller 
reflection fraction of around 1 (Hasenkopf, Sambruna & Er- 
acleous 2002). An XMM-Newton observation later revealed 
a broad Fe line, which implied an inner disc radius close 
to the innermost stable circular orbit for a maximally spin- 
ning black hole (Ballantyne & Fabian 2005) . At low energies 
the XMM-Newton spectrum showed evidence for both cold 
and warm absorption (Ballantyne 2005). 4C -1-74.26 has also 
been detected in the Smft BAT and INTEGRAL surveys 
(Tueller et al. 2007 and Bodaghee et al. 2007, respectively). 

In this paper we present the results of a 158 ks (~ 90 ks 
net exposure) Suzaku observation of 4C +74.26 performed 
in October 2007. We also carry out joint fits with the XMM- 
Newton observation performed in February 2004, and briefiy 
discuss the data from the Swift BAT and INTEGRAL cat- 
alogs. Our main goal is to try and determine the inner ra- 
dius of the accretion disc. This paper is organised as follows: 
Section [5] describes the observations and the data reduction. 
Section |3] describes the variability and Section |4] describes 
the spectral analysis. We finally discuss our results in Sec- 
tion [S] and present our conclusions in Section |S] 



2 OBSERVATIONS AND DATA REDUCTION 

4C +74.26 was observed by Suzaku between 2007 October 
28-30 for a total duration of 158 ks. Event files from version 
2.1.6.16 of the Suzaku pipeline processing were used and 
spectra were extracted using xselect. 

For each XIS, source spectra were extracted from cir- 
cular regions of 4.3 arcmin radius centred on the source 
(which was observed off-axis in the HXD nominal position) . 
Background spectra were extracted from two circular re- 
gions with the same total area as the source region, avoid- 
ing the chip corners with the calibration sources. Response 
matrices and ancillary response files were generated for each 
XIS using xiSRMFGEN version 2007-05-14 and xissimarfgen 
version 2007-09-22. The ARF generator should account for 
the hydrocarbon contamination on the optical blocking filter 
(Ishisaki et al. 2007). 



The net exposure time of all three XIS detectors 
is 92 ks and the 2-10 keV count rates are 0.753 ± 
0.003 counts s~^(XISO), 0.806+0.003 counts s~^ (XISl) and 
0.854 + 0.003 counts s"^ (XIS3). Because of the very simi- 
lar spectra of the two front-illuminated (FI) detectors (XISO 
and XIS3), their spectra were co- added for the spectral anal- 
ysis, and their response files were combined accordingly. 

For the HXD/PIN detector, a model for the non-X- 
ray background (NXB) was provided by the HXD team. 
Source and background spectra were constructed from iden- 
tical good time intervals, and the exposure time of the back- 
ground spectrum was increased by a factor of 10 (to ac- 
count for the fact that the background model was generated 
with 10 times the actual count rate in order to minimize 
the photon noise). The source spectrum was corrected for 
dead time, leaving a total net exposure time of 86 ks. The 
response file appropriate for the HXD aim point was used 
(ae_hxd_pinhxnome4_20070914.rsp) 

The total PIN count rate over the 14-30 keV energy 
range is 0.3860 + 0.0021 counts s"\ compared to 0.2743 + 
0.0005 counts s~^ for the background. Since the background 
is so much brighter than the source, the accuracy of the back- 
ground model is very important for our results. According 
to Mizuno et al. (2007), the one-sigma uncertainly of ver- 
sion 2 background models is around 3 per cent. It is often 
possible to determine the accuracy of the background model 
in a given observation more precisely by comparing it with 
the night earth spectrum. However, in the case of 4C +74.26 
the PIN field of view was never completely obscured by the 
earth, and such a comparison is not possible. Some informa- 
tion on the accuracy of the background can also be obtained 
by comparing the light curves of the source and the back- 
ground. This will be carried out in section |3] 

The background model discussed above does not include 
the contribution from the cosmic X-ray background (CXB). 
In order to account for this we simulate the spectrum of the 
CXB and add it to the spectrum of the non-X-ray back- 
ground. For the simulation of the CXB spectrum we use a 
model of the form 8 x 10~^(£/1 keV)"^'^^ exp(-E/40 keV), 
which is based on the HEAO-Al spectrum, renormalized to 
the HXD field of view. The CXB count rate in the 14-30 keV 
band is about 6 per cent of the count rate of the total back- 
ground (NXB + CXB). 

2.1 The XMM-Newton observation 

XMM-Newton observed 4C +74.26 on 2004 February 06 for 
a total duration of 34 ks. The data were reduced as described 
in Ballantyne (2005) , using the XMM-Newton Science Anal- 
ysis System version 7.1.0. We use only the data from the 
EPIC pn camera, which, after removing intervals of back- 
ground flaring, has a total good exposure time of 29 ks. 



3 VARIABILITY 
3.1 light curves 

Fig. [1] shows the 0.5-10 keV XIS light curve from the 
Suzaku observation of 4C +74.26. The light curve was ob- 
tained by adding the light curves of all three XIS detectors, 
using fully-exposed 1 ks bins. The source brightens by about 
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Figure 1. 0.5-10 keV XIS (0+1+3) light curve on a 1 ks time- 
scale, using only fully-exposed bins. 
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Figure 2. 14—30 koV PIN light curves on a 20 ks time-scale. 
Panels from top to bottom show the total light curve, the light 
curve of the NXB model, and the source light curve (where both 
the NXB model and the constant contribution from the CXB have 
been subtracted) . As a reference the bottom panel also shows the 
XIS light curve (dotted line), rescaled to the mean PIN count 
rate. The background-subtracted light curve is clearly correlated 
with the background model, suggesting that the background has 
been underestimated. 



20 per cent during the observation. We also note that the 2- 
10 keV flux has gradually increased by about a factor of two 
since the source was observed by ASCA in 1996, as summa- 
rized in Table [1] The flux in the soft band in the ASCA ob- 
servation was also found to be two times higher than in a 
ROSAT observation in 1993 (when the source was in the 
field of view of an observation targeted at the cataclysmic 
variable VW Cep, Brinkmann et al. 1998). The 2-10 keV 
fiux in the Suzaku observation is 3.24 x 10~^^ erg cm~'^ s~^, 
corresponding to a luminosity of 8.80 x 10^** erg s~^. 

In Fig [2] we show 14-30 keV PIN light curves on a 20 ks 
time-scale. Specifically, we compare the total light curve, 
the light curve of the NXB model, and the light curve of the 
source. The source light curve was obtained by subtracting 
both the NXB model and the constant contribution from 
the CXB. As a reference we also show the XIS light curve, 
rescaled to the mean PIN count rate. The light curves of 
the background model and the source are clearly correlated 



Obs. date 


Mission 


Z-IV KeV nux 


Reference 






(10""'^"'^ erg cm~^ 




2007-10-28 


Suzaku 


3.24 


(1) 


2004-02-06 


XMM-Newton 


2.43 


(2) 


1999-05-17 


BepppoSAX 


1.41 


(3) 


1996-09-09 


ASCA 


1.69 


(4) 



Table 1. 2-10 keV fluxes of 4C +74.26. References: (1) This pa- 
per, (2) Ballantyne & Fabian (2005), (3) Hasenkopf et al. (2002), 
(4) Sambruna et al. (1999). 
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Figure 3. hardness ratios as a function of count rate for XIS and 
3 (black and blue crosses, upper panel) and XISl (lower panel). 
The hardness ratio is defined as HR=(H-S)/(H+S), where S is 
the count rate in the 0.5 — 2 keV band and H is the count rate in 
the 3 — 10 keV band. The hardness ratios are calculated based on 
~ 6 ks spectra. The thin crosses represent un-binned data points 
and the thick ones show the mean of data points binned along the 
X-axis. Each bin contains four original data points, apart from the 
last bin which contains three. For all three XIS we see that the 
spectrum tends to harden as the count rate increases. 



(the cross-correlation function peaks at a value of 0.82 at 
time lag), indicating that the background model has been 
underestimated. 
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Figure 4. Rms spectrum of the entire observation calculated on 
the orbital time scale (~ 6 ks). The spectrum was calculated as 
an average of the rms spectra of the three XIS detectors. 



3.2 Hardness ratios 

In order to characterise the spectral variabihty of the source 
we start by calculating hardness ratios, HR = H — S/H + S, 
where we take H to be the count rate in the 3-10 keV band 
and S to be the count rate in the 0.5-2 keV band. The 
hardness ratios for all three XIS, calculated in orbital length 
bins of 5760 s, are plotted as a function of the H + S count 
rate in Fig. [31 In order to reduce the scatter in the plot we 
binned consecutive points along the x-axis so that each bin 
comprises four original data points (apart from the last bin 
which contains three data points). Un-binned and binned 
data are shown as thin and thick crosses, respectively. 

For all three XIS we see that the spectrum tends to 
harden as the count rate increases, although within a fairly 
large spread. In agreement with this, we find that a positive 
correlation is a better fit than a constant for all three data 
sets. According to an F-test, the improvement in the fit for 
un-binned (binned) data is significant at 91 (93) per cent 
(XISO), 99 (99) per cent (XISl) and 98 (87) per cent (XIS3). 
The fact that the positive correlation is highly significant 
only in XISl is presumably due to the higher effective area 
at low energies in this detector, which reduces the error on 
the count rate in the soft band. 



3.3 The rms spectrum 

We next consider the root-mean-squared (rms) spectrum, 
which, for a given time-scale, shows the fractional variabil- 
ity as a function of energy. The techniques for calculating 
rms spectra are described in e.g. Edelson et al. (2002) and 
Vaughan et al. (2003). Fig. |4] shows the rms spectrum of 
4C -1-74.26, calculated as an average of the rms spectra of 
the three XIS detectors, on the orbital time-scale of 5760 s. 
The errors from the Poisson noise were calculated following 
Vaughan et al. (2003). The rms variability clearly increases 
with energy up to about 4 keV, at which point it flattens 
out and shows some evidence of decreasing. 




Energy 

Figure 5. FI XIS -I- PIN spectrum as a ratio to a power law (mod- 
ified by Galactic absorption), fitted over the 2—4 and 7-10 keV 
energy ranges. The red and green data points show the PIN data 
assuming a background model of ±3 per cent (the one-sigma ac- 
curacy of v2 background models). The study of the PIN light 
curves in section 13.11 showed that the background has probably 
been underestimated, suggesting that the real value should be 
closer to the red data points. The low-energy data clearly reveal 
the presence of absorption in excess of the Galactic value. 



4 SPECTRAL ANALYSIS 

In the spectral analysis of the Suzaku data we will concen- 
trate on the 2-10 keV (XIS) and 14-30 keV (PIN) energy 
ranges. We find that the agreement between the FI XIS is 
very good in the 2-10 keV range, and we therefore add their 
spectra, as previously mentioned in section [2] and recom- 
mended in the Suzaku Data Reduction Guide. The spectrum 
of XISl, on the other hand, differs significantly from the FI 
XIS spectra in the Fe Ka region (the line energy and width 
are inconsistent with the FI XIS values), and we do not in- 
clude it in our spectral analysis. Because of inconsistencies 
between the XIS detectors at low energies (presumably due 
to the modelling of the contamination), we do not extend 
the spectral analysis below 2 keV at this time. We note, 
however, that the low-energy spectrum is characterized by 
absorption in excess of the Galactic value (see Fig. [SJ, in 
agreement with previous observations. The absorption is not 
strong enough to affect the spectrum above 2 keV. 

Below we will first consider the 2-10 keV FI XIS spec- 
trum and then include the high-energy PIN data. We will 
also investigate the high-energy properties of the source us- 
ing the spectra from Swift and INTEGRAL. As a last step 
we will perform joint fits to the Suzaku and XMM-Newton 2- 
10 keV spectra. All spectral fits were performed using xspbc 
version 11.3.2aj. Errors on model parameters are quoted at 
the 90 per cent confidence level and energies of spectral fea- 
tures are quoted for the rest frame of the source. All fits 
include Galactic absorption fixed at Nh = 1.19 x 10^^ cm~^ 
(Dickey & Lockman 1990). 



4.1 The 2-10 keV FI XIS spectrum 

In order to look for the presence of a broad Fe line we start 
by fitting the spectrum with a power law, excluding the 4- 
7 keV energy band, where the line emission should dominate. 
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Figure 6. Confidence contours for rin and i in the model consist- 
ing of a power law and a Laor line. The black (solid), red (dashed) 
and green (dotted) lines represent 1, 2 and 3 sigma confidence, 
respectively. The best fit is shown as a black dot. 



Fig. E] shows the data as a ratio to the power-law fit (which 
has r = 1.81 ± 0.01), clearly revealing that a broadened line 
is present. The fact that the line is broadened is confirmed 
when we add a Gaussian line to the power-law model. If 
we let the width of the line be a free parameter we find an 
excellent fit of = 1510 for 1509 de grees of freedom (d.o.f.), 
compared to x^/d.ol. = 1546/1510 if we force the line to be 
narrow {a = 1 eV). The broad line has a = 0.24j;[;:^? keV 
and EW = 85tig eV. The results of these fits as well as all 
other spectral fits discussed in this section are presented in 
Table [II 

Broad Fe lines are often found to have a narrow com- 
ponent originating from distant matter. To check for such a 
component in 4C -1-74.26 we add a narrow line at 6.4 keV 
to the model with the broad Gaussian. This improves the 
fit by only Ax^ ~ —2 for 2 fewer degrees of freedom, and 
the narrow line is found to have a very low equivalent width 
of lOtJo We thus conclude that a narrow component to 
the Fe line is not required. 

We next replace the broad Gaussian with a Laor line 
(Laor 1991), which models the Fe line emission from an 
accretion disc around a Kerr black hole. We fix the en- 
ergy of the line at 6.4 keV, the outer radius of the disc at 
Tout = 400 rg (where rg = GM/c? is the gravitational ra- 
dius), and the emissivity index of the disc at g = 3 (the 
emissivity follows the form e oc r"', where r is the ra- 
dius of emission, and g = 3 is expected from gravitational 
energy release). This leaves the inner radius of emission 
(rin), the inclination of the disc (i) and the normalization 
of the line as free parameters. We find an excellent fit of 
xVd.o.f. = 1514/1509 with = 53^28^ rg and i = 20t®°. 

The relatively low inclination obtained in this fit implies 
that the source is extremely large (of the order of 3 Mpc). 
This is still smaller than the largest radio galaxy observed to 
date (3C 236 is ~ 4.5 Mpc, Saripalh & Mack 2007) but larger 
than the majority of the known radio galaxies. Assuming 
that the size of 4C -1-74.26 is around 2 Mpc, which is more 
typical for a giant radio galaxy, the inclination should be 
around 33° . If we fix the inclination at this value we obtain 
a fit which is only slightly worse (xVd-o.f. = 1518/1510), 
but which has a significantly larger rin of 177lyy^ rg. The 
relationship between the inclination and rin can be seen in 



Fig.[6l which shows the confidence contours for these param- 
eters in the fit where the inclination was left free. We note 
that the three-sigma confidence contour remains below 50° , 
in agreement with the upper limit on the inclination of 49° 
set by Doppler boosting arguments (Pearson et al. 1992). 

The relatively large values obtained for rin in the fits 
described above suggest that the inner disc in 4C +74.26 
is either not present or for some reason unable to produce 
strong fluorescent Fe emission. In order to check if a fit where 
the Fe emission is produced in the inner regions of the disc 
can be obtained, we refit the data with rin fixed at 1.24 rg 
(the innermost stable circular orbit of a Kerr black hole). 
We find that a narrow 6.4 keV line has to be added to the 
model in order to obtain an acceptable fit for this scenario. 
The quality of the fit is worse than for the case of a free 
rin (Ax^ = +15 for one less d.o.f) but we note that the 
combination of the small rin and the narrow line results in 
a higher inclination of 40^6°. 

In order to also model the reflection spectrum associ- 
ated with the Fe line, we next replace the Laor line with 
the REFLiON model by Ross & Fabian (2005), which mod- 
els the emission from a constant-density accretion disc. The 
parameters of the model are the Fe abundance, the ioniza- 
tion parameter 5, the photon index of the incident power 
law and the normalization. In order to account for the rela- 
tivistic effects in the vicinity of the black hole, the reflection 
model is convolved with the relativistic blurring kernel kd- 
BLUR, which is derived from the code by Laor (1991) (the 
model parameters are the inner and outer radius of the disc, 
the inclination and the emissivity index). We leave all pa- 
rameters free to vary apart from the outer radius of the disc 
and the emissivity index, which we fix at rout = 400 rg and 
q ~ S, respectively. 

We find a best fit of xVd-o.f. = 1514/1506 with nn = 
68^25^ rg and i = 2llg*°, which is very similar to our best- 
fitting Laor model. As in the case of the Laor model, the fit 
only becomes slightly worse if the inclination is fixed at a 
larger value. We are also still able to obtain an acceptable 
fit with rin fixed at 1.24 rg if a narrow line at 6.4 keV is 
included. As before, this fit is somewhat worse (Ax^ = +13 
for one less d.o.f) and has a higher inclination (42j;5°). In 
both reflection fits the ionization parameter of the disk is 
30 erg cm s~^, i.e nearly neutral, and the Fe abundance is 
about 2 times solar. Both models are plotted in Fig. [7] to- 
gether with the flt residuals. It is clear that the first model 
(with the larger rin and lower inclination) is a better fit in 
the Fe line region. We also note the presence of a narrow 
feature around 7.2-7.4 keV (6.5-6.7 keV in the observed 
frame). We will discuss the possible origin of this feature in 
section 14.1.11 below. 

In addition to the interpretation that the inner disc is 
missing, there are several possible scenarios that could ex- 
plain the large values obtained for rin in our best-fitting 
models. For example, extreme ionization of the inner disc 
would lead to a nearly featureless reflection spectrum that 
would be hard to disentangle from the power law contin- 
uum. In fact, adding an inner reflector with ^ fixed at 
10* erg cm s~^ (the highest allowed value) to our best-fitting 
reflection model results in an acceptable flt, although the 
quality of the data is not good enough to constrain the fit 
parameters. 

Another possible reason for the large value of rin is that 
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Table 2. Fits to the 2—10 keV Suzaku FI XIS spectrum. The model components are pi = power law, nl = narrow Gaussian emission 
line ((T fixed at 1 eV), G = Gaussian emission line, Laor = Laor model, blr = KDBLUR (relativistic blurring kernel), ref = REFLION 
(reflection model). F is the photon index of the power law, Eq/l is the energy of the broad Gaussian or Laor line (in keV), ctq is the 
width of the Gaussian line (in keV), EWg/l the equivalent width of the broad Gaussian or Laor line (in eV), Eni is the energy of the 
narrow line (in keV), EWni is the equivalent width of the narrow line (in eV), r-^^ is the inner radius of the disc (in gravitational radii), 
i is the inclination of the disc to the line of sight (in degrees), Fe is the iron abundance of the disc (relative to solar), ^ is the ionization 
parameter of the disc (in erg cm~-^ s~^). Superscript f indicates that the parameter was frozen. The following parameters were always 
fixed: a of the narrow line (1 eV), the outer radius of the disc (400 gravitational radii), the emissivity index of the disc (3). All fits also 
included absorption fixed at the Galactic value. 
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Figure 7. Reflection fits to the 2—10 keV Suzaku FI XIS spectrum. 
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The top panel shows the best-fitting models together with their 
components, and the bottom panel shows the residuals to the fit. The model in the left panel consists of a power law and reflection from 
a disc truncated around 70 rg. The model shown in the right panel consists of a power law, reflection from a disc that extends all the 
way in to 1.24 rg, and a narrow Fe Ka line. Black (solid) lines represent the total model, red (dashed) lines the power law, blue (dotted) 
lines the reflection model, and the green (dashed-dotted) line the narrow line. 
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the illumination of the disc results in an emissivity index 
that deviates from our assumed value of g = 3. The quality 
of the data does not allow us to constrain this parameter if 
all other parameters are left free. However, when fixing the 
other reflection parameters at their time-averaged values, we 
find that the model with rin = 1.24 rg improves significantly 
with a flatter emissivity of g ~ 1. Such a flat emissivity 
would be expected (at least in an average sense) if the disc 
is illuminated from a great height (> 30 rg) in a lamp post 
geometry (see also Vaughan et al. 2004), or if the power-law 
emission originates in the base of a jet, which is beamed 
away from the inner disc. 

A contribution from a jet in 4C 4-74.26 has previously 
been suggested for the ASCA observation of the source 
(Sambruna et al. 1999). In order to model this scenario we 
add a second, hard power law to the reflection models de- 
scribed above. The results of this are inconclusive though, 
as we find that we cannot obtain any meaningful constraints 
on the slope and normalization of the second power law. 

4-.1.1 The narrow feature around 7.4 keV 

As mentioned above, the residuals in Fig. [7] reveal a narrow 
feature around 7.2-7.4 keV (6.5-6.7 keV in the observed 
frame). Depending on the underlying model, this feature 
can be modelled either as a 7.40 keV emission line or as a 
7.17 keV absorption line. 

For the model consisting of a power law and a broad 
Gaussian, the best fit is obtained with a narrow emis- 
sion line at 7.40 keV (the fit improves by A.x^ = —8 for 
2 fewer degrees of freedom when the line is added, and 
EW= 16 ± 11 eV). The same is true for the Laor and re- 
flection models in table [2] that have rin as a free parame- 
ter. None of the best-fitting parameter values change signifi- 
cantly when the line is added to these models. An absorption 
line at 7.17 keV results in a very small improvement of the 
fit for these models, and we especially note that the addition 
of the line does not significantly affect the inclination. This 
is also true if the absorption line is identified with Fe xxvi 
at 6.97 keV and we add a Fe xxv line (rest energy 6.7 keV) 
with the same blueshift. 

For the models where rin = 1.24 rg and a narrow 
6.4 keV line is included, we find that the narrow feature 
can be equally well fitted with a 7.17 keV absorption line 
as with a 7.40 keV emission line. Specifically, in the case 
of the model which consists of a power law, a Laor line 
and a 6.4 keV line, we find that the absorption line im- 
proves the fit by Ax^ = — 13 for 2 fewer degrees of free- 
dom and has EW= 22^]^^ eV, while the emission line re- 
sults in Ax^ = — 14 for 2 fewer degrees of freedom and has 
EW= 16ljj eV. No model that includes both an absorption 
and an emission line can be found. Although the absorption 
line is an equally good fit as the emission line in this case, 
we note that the inclination goes up to about 55° when the 
absorption line is added. This is higher than the upper limit 
of 49° placed by Doppler boosting arguments (Pearson et 
al. 1992). 

We have ruled out an instrumental origin for the 
lines, which suggests that they (if real) are due to highly 
blueshifted Fe emission. For example, identification of the 
7.17 keV absorption line with Fe xxv at 6.70 keV would im- 
ply an outflow velocity of about 20 000 km s~^. This extreme 



velocity together with the relatively low signiflcance of the 
lines suggests that they are probably not real. However, it 
is interesting to note that Robinson et al. (1999) detected a 
redshifted Ha line in polarised light in this source, indicating 
that an outflow with a velocity greater than 5 000 km s~^ 
is present. 



4.2 The high-energy PIN spectrum 

4C -1-74.26 is clearly detected in the PIN up to about 30 keV 
and we stress that this detection is robust against uncertain- 
ties in the background model. When fltting the data we use 
the model for the non- X-ray background provided by the 
HXD-team, and we account for the contribution from the 
CXB as described in section (2] The cross normalization of 
the PIN with respect to XISO has been reported to be about 
1.14 for the HXD nominal position and v2 data (Ishida, 
Suzuki & Someya 2007). After adjusting this value for the 
difference in normalization between XISO and our added FI 
XIS spectrum, we adopt a PIN/FI XIS cross normalization 
of 1.12. 

Fig. [5] shows the PIN -f FI XIS spectrum as a ratio to a 
power law, fitted over the 2-4 and 7-10 keV energy ranges. 
In order to illustrate the uncertainties in the PIN data we 
show the spectra obtained by increasing and decreasing the 
background model by 3 per cent (the one-sigma uncertainty 
of v2 background models). We note that the study of the 
PIN light curves in section ISTTI showed that the background 
has probably been underestimated, suggesting that the real 
value is closer to the lower (red) data points. This also agrees 
with the predictions of the reflection models that were found 
to fit the 2-10 keV data. For the model where rin is a free pa- 
rameter we get a best 2-30 keV fit of xVd-o.f. = 1561/1549 
if the background is increased by 2 per cent. Similarly, for 
the model with rin = 1.24 rg, we get a best chi-squared of 
X'^/d.o.f. — 1573/1548 if the background is increased by 1-2 
per cent. 

The relative contribution from the reflection component 
is usually measured by the reflection fraction, which is the 
ratio of the total reflected emission to the power-law emis- 
sion. For an isotropic source illuminating an inflnite plane 
of gas the reflection fraction is 1. We estimate the reflection 
fraction of our models by measuring the the unabsorbed 
fluxes in the two components within xspec. For the two 
models presented above the reflection fractions are 0.3 and 
0.5, respectively, assuming that the power law emission ex- 
tends all the way up to the upper boundary for the reflion 
model. On the other hand, assuming an upper cutoff for the 
power law of 150 keV (as previously found for BepposSAX 
data, Grandi et al. 2006) we find slightly higher reflection 
fractions of 0.4 and 0.7. In either case, it is clear that the 
reflection fraction is less than 1, implying e.g. that the il- 
lumination is anisotropic, that the disc is truncated or that 
an unreflected component (such as a jet) contributes to the 
spectrum. 

4.3 Comparison with the Swfit and INTEGRAL 
data 

4C -1-74.26 is included in both the Swift BAT 9-month AGN 
catalog (Tueller et al. 2007) and the INTEGRAL general 
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Figure 8. The Suzaku FI XIS spectrum (in black), together with 
the Swift (red) and INTEGRAL (blue) spectra, shown as a ratio 
to the best-fitting 2-10 keV reflection model. The normalizations 
of the Swift and INTEGRAL data were allowed to vary freely 
with respect to the XIS spectrum. 
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Figure 9. Comparison of the XMM-Newton EPIC pn (black) 
and Suzaku FI XIS (red) spectra. The plot shows the spectra as a 
ratio to a power law fitted to the XMM-Newton data, excluding 
the 4—7 keV energy band. 



reference catalog (Bodhagee et al. 2007). In order to fur- 
ther investigate the high-energy properties of 4C -1-74.26 we 
briefly investigate the spectra from these catalogs in this sec- 
tion. The data from the BAT catalog are publicly availabl^B 
and the INTEGRAL spectrum was provided by the INTE- 
GRAL AGN team (the INTEGRAL data is also analysed in 
Molina et al. 2008, in prep). 

A simple power-law fit to the data yields T = 2.1^03 
for Swtft (14-195 keV) and F = 2.2l°;^ for INTEGRAL 
(20-100 keV). This should be compared with T = 1.8 ± 0.2 
in the 14-30 keV range for the Suzaku PIN detector. No 
clear evidence for a cutoff of the power law is seen in any 
of the spectra. Based on the power-law fits the inferred 
14-30 keV fiuxes are ^^14-30 = 1-4 x 10"" erg cm"^ s"^ 
(Swift) Fi4_30 = 2.6 X 10"" erg cm'^ s"^ (INTEGRAL) 
and F14-30 = 2.7 x 10~" erg cm"^ s"^ (Suzaku PIN). As 
the observations are not simultaneous and the source ap- 
pears to be getting brighter (as shown in Table [T|, it is not 
surprising that the photon indices and fluxes differ some- 
what. 

In order to check if the Swift and INTEGRAL data 
are consistent with a reflection model, we flt them with 
our best-fitting 2-10 keV refiection model, leaving only the 
cross-normalizations between the detectors as free param- 
eter. As shown in Fig. [S] the shape of the high-energy 
spectra can be well reproduced by the reflection model 
(x^/dof = 1527/1524 over the entire 2-195 keV range). The 
cross-normalizations with respect to the FI XIS are found to 
be 0.6 for Swift and 0.9 for INTEGRAL. We also note that 
the high-energy data can be equally well described with the 
reflection model with rin = 1.24 Vg. 

4.4 Fits including the previous XMM-Newton data 

Fig. |9] shows a comparison of the Suzaku FI XIS spec- 
trum and the EPIC pn spectrum from the previous XMM- 
Newton observation of 4C -1-74.26 (Ballantyne & Fabian 

^ http:/ /swift. gsfc.nasa.gov/docs/swift/results/bs9mon/ 



2005). The spectra are shown as a ratio to a power- law fit- 
ted to the XMM-Newton data, excluding the 4-7 keV energy 
range. It is clear from the plot that the continuum is some- 
what steeper in the Suzaku observation (AF ~ 0.1), while 
the Fe line profiles are very similar in both spectra. 

In order to to try and find a spectral model for both 
observations we performed simultaneous fits to the spec- 
tra, using the two different reflection fits obtained for the 
Suzaku data as a starting point. Experimentation showed 
that very good joint fits could be obtained for both models 
if all parameters apart from the slope and the normalization 
of the power law were tied between the two spectra. The 
results of the fits are presented in Table |3] We see that the 
model where rin is a free parameter is still a somewhat better 
fit than the model with n,, = 1.24 rg (xVd-o.f. = 2410/2538 
compared to 2431/2537). In both models the photon index 
of the power law is 0.1 steeper for Suzaku than for XMM- 
Newton . The only tied parameters that differ significantly 
from the values obtained when fitting only the Suzaku data 
are rin, i and Fe of the first model. The values of these pa- 
rameters have all increased (especially note that the inclina- 
tion now has a more reasonable value of 3llij°) but they are 
all still consistent (within the errors) with the Suzaku fits. 
As in the case of the Suzaku fits, we also find that both mod- 
els agree well with the high-energy PIN data, and that the 
best agreement is found if we assume that the background 
model has been underestimated by 2-3 per cent. 

We stress that the reflection normalization is tied in 
these flts and that letting it vary independently for the two 
spectra does not improve the quality of the fits. The same 
is true for the normalization of the narrow 6.4 keV line. We 
thus conclude that the difference between the Suzaku and 
XMM-Newton observations can be well described in terms 
of a power law that has increased in normalization and also 
become slightly steeper. 

As a final comment we note that there is evidence 
for a narrow 6.2 keV line in the XMM-Newton spectrum 
(EW ~20 eV, as previously noted by Ballantyne & Fabian 
2005) but not in the Suzaku spectrum. The origin of this line 
is unclear but Ballantyne & Fabian (2005) suggested that it 
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pi + blr*rcf 

pi + blr*rcf + 111 


1 83+0-02 


„ +0.01 
^•'^-0.02 

1 76+0'05 

^- '"-0.05 


j^gj^+239 

1.24' 


311- 
401^ 


3.8l|;2 


30+21 
30+0 


6 4+0.04 
"■^-0.02 


281^ 


36+10 


2410/2538 
2431/2537 



Table 3. Combined 2—10 keV reflection fits to the Suzaku and XMM-Newton spectra. The model components are pi = power law, blr = 
KDBLUR (relativistic blurring), ref = REFLION (reflection model), nl = narrow Gaussian emission line (cr fixed at 1 eV). F is the photon 
index of the power law, rin is the inner radius of the disc (in gravitational radii), i is the inclination of the disc to the line of sight (in 
degrees), Fe is the iron abundance of the disc (relative to solar), ^ is the ionization parameter of the disc (in erg cm~2 s~l), E^i is the 
energy of the narrow line (in keV), EWni is the equivalent width of the narrow line (in eV). Subscripts s and x denote values for the 
Suzaku and XMM-Newton spectra, respectively. All parameters were tied between the two spectra apart from the photon index and the 
normalization of the power law. Superscript f indicates that the parameter was frozen. The following parameters were fixed in all the 
fits: the outer radius of the disc (400 gravitational radii), the emissivity index of the disc (3), a of the narrow line (1 eV). All fits also 
included absorption fixed at the Galactic value. 



and the 6.4 keV line could form the red and blue horns of a 
diskline originating in the outer disc. 



5 DISCUSSION 

5.1 Spectral variability 

Due to the long duration of the Suzaku observation, we have 
for the first time detected variability in the light curve of 
4C +74.26. We find that the source brightens by about 20 
per cent during the observation, implying a doubling time of 
~ 1 Ms. There is evidence for a positive correlation between 
the hardness ratio and the count rate, although this is only 
significant at the 99 per cent level in the case of XISl. We 
also find that the rms variability increases with energy up 
to about 4 keV, at which point it flattens out and shows 
evidence of decreasing. This spectral variability is rather 
different from that observed in typical radio-quiet Seyfert 
galaxies, which generally exhibit spectral softening with in- 
creasing count rate and rms variability decreasing with en- 
ergy above ~ 1 keV (e.g. Papadakis et al. 2002; Markowitz 

6 Edelson 2004). 

One possible explanation for the spectral variability ob- 
served in 4C -1-74.26 is that a jet contributes to the spectrum. 
This is possible even if the inclination of the jet to the line 
of sight is fairly large, as the base of the jet may not be well 
coUimated. In this picture, our measured photon index of 
~ 1.8 would be due the combined contributions from a hard 
component associated with the jet, and a softer component 
associated with Comptonized emission from the corona, as 
previously suggested for an ASCA observation of the source 
(Sambruna et al. 1999). The fact that the source is harder 
when it is brighter would then indicate that the relative 
contribution from the jet is higher when the flux increases, 
assuming that the intrinsic slopes of the two power laws 
do not vary significantly. Such a model could in principle 
be tested by investigating the difference spectrum between 
high- and low-fiux spectra. However, due to the limited am- 
plitude of the variability in 4C -1-74.26 we find that the dif- 
ference spectrum is not of sufficient quality to discriminate 
between models. 



5.2 Spectral modelling 

In our spectral analysis of the Suzaku data of 4C 4-74.26 we 
clearly detect a broad Fe line, as previously seen by XMM- 
Newton (Ballantyne & Fabian 2005). Our best fits to the 
Fe line profile suggest that the line originates outside about 
50 Tg and that the inclination to the line of sight is around 
20° . This inclination implies that the source is around 3 Mpc 
across, which is extremely large even for a giant radio galaxy. 
Although it is possible that the source actually is this large, 
we note that an almost equally good fit can be obtained if 
the inclination is fixed at a larger value. For this scenario the 
inner radius of emission is outside 100 Tg. We also managed 
to find a good fit with a model in which rin was fixed at 
1.24 Tg. This model has the advantage of a higher inclination 
of around 40°, but is a worse fit in the Fe line region and 
also relies on the presence of a narrow Fe Kq line, which is 
not required in any other model. Overall it seems like the 
first model is a better approximation of the real conditions. 

The perhaps simplest explanation for these findings is 
that the inner accretion disc is missing, as suggested for 
many other radio-loud AGN (e.g. Eracleous et al. 2000). 
It is also possible, however, that the inner disc is highly 
ionized or not strongly illuminated. The latter explanation 
seems more likely, as the former model would require a very 
sharp transition from an extremely ionized to a nearly neu- 
tral disc. A plausible explanation for why the inner disc is 
not strongly illuminated is that the corona, where the power 
law originates, is located high above the disc. 

The coronae of AGN are widely believed to be powered 
by magnetic fields in their accretion discs. As magnetic fields 
are thought to be important for jet formation (e.g. Blandford 
& Payne 1982), it is likely that the properties and structure 
of the magnetic fields in radio-quiet and radio-loud objects 
are different. The strong reflection features often seen in 
radio-quiet AGN indicate that their inner discs are strongly 
illuminated, implying that the magnetic energy is dissipated 
close to the disc. It is conceivable that the magnetic field 
structure in radio-loud objects leads to the energy being dis- 
sipated at larger heights, possibly in a region moving away 
from the disc, thus resulting in less illumination of the inner 
disc and weaker refiection features. 

It is also possible that emission from the jet itself con- 
tributes to the spectrum, as tentatively suggested by the 
variability properties described above. We explicitly tested 
for the presence of a jet in the Suzaku spectrum by adding a 
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second, hard power law to the power law + reflection model. 
The results were inconclusive, however, as we found that we 
could not constrain the properties of the second power law. A 
much longer observation would be needed to put constraints 
on the possible jet component in 4C +74.26. 

The difference between the 2007 Suzaku and 2004 
XMM-Newton spectra can be well modelled with a power 
law that has increased in normalization and become slightly 
steeper. The reflected emission is consistent with staying 
constant between the observations, and the best-fitting disc 
parameters in the joint fits are similar to those found when 
only fitting the Suzaku data. A possible explanation for this 
behaviour is that the increase in flux between the two obser- 
vation is mainly due to the (unreflected) emission from the 
jet increasing, in combination with the slope of the coronal 
and/or jet power laws steepening somewhat. 



6 SUMMARY AND CONCLUSIONS 

We have analysed a 92 ks exposure Suzaku observation of 
the radio quasar 4C 4-74.26. Our main results can be sum- 
marized Eis follows: 

• The 2-10 keV flux measured by Suzaku is 
3.24 X 10~^^ erg cm~^ s~^, which is the highest flux 
observed to date. The source has gradually brightened by 
about a factor of two since it was first observed by ASCA in 
1996. 

• Due to the long observation, we have for the first time 
detected variability in the light curve of 4C -1-74.26. We find 
that the count rate increases by about 20 per cent during 
the observation. 

• There is evidence for spectral hardening as the 
count rate of the source increases. We also find that 
rms variability increases with energy up to about 4 keV, 
where it flattens out and shows some evidence of decreasing. 

• We clearly detect a broadened Fe Ka line. When 
fitted with a simple Gaussian, it has EW = 85 eV and 
a = 0.24 keV. 

• The entire 2-30 keV energy range can be well fitted 
with a power law (F ~ 1.8) and near-neutral reflection 
from an accretion disc. No reflected emission inside 
about 50 Tg is required by the Fe line proflle, although 
a model with rin = 1.24 rg provides a good flt if a 
narrow Fe Ka line is included. The reflection fraction is 
found to be in the range ^ 0.3-0.7, depending on the model. 

• The high-energy Swift BAT and INTEGRAL spectra 
show a somewhat steeper continuum slope F ~ 2.1 — 2.2) 
but can nevertheless be well described with the power-law 
-I- reflection models described above. 

• The difference between the Suzaku and XMM- 
Newton 2-10 keV spectra can be well explained with a power 
law that has increased in normalization and become slightly 
steeper. 



The spectral variability properties described above can 
be explained in a model where the relative contribution from 
a hard jet increases with increasing flux. However, because of 
the limited variability amplitude and spectral quality, we are 
not able to constrain the possible contribution from a jet in 
4C +74.26. The best-fitting reflection model implies that the 
inner disc is either truncated at large radii or not strongly 
illuminated. We suggest that the latter scenario may be due 
to the power-law source being located high above the disc, 
e.g. as a result of the magnetic field structure. 
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